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Abstract: Low-coherence spectroscopy (LCS) offers the valuable 
possibility to measure quantitative and wavelength resolved optical property 
spectra within a tissue volume of choice that is controllable both in size and 
in depth. Until now, only time domain detection was investigated for LCS 
(tdLCS), but spectral domain detection offers a theoretical speed/sensitivity 
advantage over tdLCS. In this article, we introduce a method for spectral 
domain detection in LCS (sdLCS), with optimal sensitivity as a function of 
measurement depth. We validate our method computationally in a 
simulation and experimentally on a phantom with known optical properties. 
The attenuation, absorption and scattering coefficient spectra from the 
phantom that were measured by sdLCS agree well with the expected optical 
properties and the measured optical properties by tdLCS. 
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1. Introduction 

Low-coherence spectroscopy (LCS) is a method for simultaneous measurements of local 
attenuation, absorption, scattering and backscattering coefficient spectra in turbid media 
within the visible wavelength range (480 to 700 nm) [1-3]. LCS distinguishes itself from 
other quantitative spectroscopic techniques by controlling the photon path length of the 
detected light, which offers the unique possibility to control the measurement volume both in 
size and in depth — e.g. for the in vivo measurement of hemoglobin concentrations within 
distinct regions of the human dermal microcirculation [3]. Our previous studies on LCS 
involved a time-domain detection setup (tdLCS). However, from optical coherence 
tomography (OCT) studies, we know that spectral domain detection allows for more sensitive 
measurements, or higher acquisition speed at unaltered sensitivity [4,5]. Since clinical 
applications of LCS require both fast and sensitive measurements, we investigate the 
possibility for spectral domain detection in LCS (sdLCS), focusing primarily on validating a 
new method of analysis that was developed for this purpose. 

Signal acquisition in sdLCS is comparable to that in spectrometer-based, spectral domain 
OCT (sdOCT), with the major difference that in sdLCS optical property spectra are derived 
from the spectral content of the signal, and in sdOCT an image is reconstructed from the 
backscattered intensity. To obtain optical property spectra within confined tissue volumes, 
sdLCS requires both high spatial and high spectral resolution. The spatial resolution for the 
acquisition of one backscattered spectrum in sdLCS and sdOCT is limited by the axial 
measurement range of the spectrograph, which depends on the spectral resolution of the 
spectrograph. Since the standard methods for signal acquisition in sdOCT commonly result in 
only one backscattered spectrum within the entire axial measurement range of the 
spectrograph (typically 1 to 2 mm) [4,5], the spatial resolution of these methods is not 
sufficient for sdLCS (typically 22 |im) [1-3]. Spectroscopic sdOCT methods solve this 
problem typically by windowed Fourier transforms on the measured spectrum, which result in 
multiple depth resolved spectra [6-8]. Due to the sensitivity decrease with depth inherent to 
these methods and the absence of focus tracking, these methods may not achieve optimal 
sensitivity in depth. In addition, precise refractive index-dependent corrections for the 
unwanted signal attenuation due to these effects are needed, in order to achieve quantitative 
measurements of the tissue attenuation. Therefore, we developed a method for sdLCS that 
applies focus tracking and measures with maximum spectrograph sensitivity (around zero 
delay) at every depth. 

Our approach involves acquisition of local back scattered spectra within an axial 
measurement range of only 22 |im. Translation of this axial measurement range through the 
sample results in depth-resolved backscattered spectra acquired around zero delay with high 
spatial and sufficiently high spectral resolution (6 nm) while focus tracking is applied. The 
main challenge in our approach to sdLCS is to isolate the sample's backscattered spectrum 
from the detected spectrum, in order to retrieve the spectral absorption and backscattering 
features from the sample only. For this purpose, the unwanted DC and modulation 
components in the detected spectrum can be efficiently removed by introducing a phase 
modulation in the reference arm [9,10], for which we use an oscillating reference mirror. 

Since the primary aim of this study is to validate whether we can use this new approach 
for the determination of local absorption and scattering coefficient spectra, we will provide 
proof of principle with an on-the-shelf spectrograph (USB4000, Ocean Optics, USA). Only 
part of the feasible speed advantage will be demonstrated, since the quantum efficiency and 
acquisition speed of this spectrograph are not optimized for sdLCS. We will validate our 
algorithm with a simulation of an sdLCS absorption measurement, while simulating the 
USB4000 as the detecting spectrograph. Subsequently, we experimentally demonstrate the 
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ability of sdLCS to measure the attenuation, absorption and scattering coefficient spectra from 
a polystyrene-dye phantom, and we compare our results to those of a tdLCS measurement. 



AR=v R /(2f R ) 
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Fig. 1. Schematic illustration of setup and parameters in sdLCS/tdLCS. Xq: center wavelength, 
A,bw: spectral bandwidth, x R : reference arm length, x s : sample arm length, AOPL: path length 
difference, AR: reference mirror scanning range, v R : reference mirror scanning velocity, f R : 
reference mirror scanning frequency, I: geometrical path length in sample, Ai R : scanning 
window in sample induced by AR, n: refractive index, |u t : attenuation coefficient, i D : detector 
current, t: integration time, N p : number of pixels, 3X: pixel width, AX: wavelength resolution. 



2. Theory 

2.1. Background 

For sdLCS, we use a setup that is based on a Michelson interferometer and a broad band light 
source (Fig. 1). A spectrograph detects the spectrum i D , which (for a single reflector) is 
modulated by the optical path length difference between the sample arm and the reference arm 
AOPL = 2(x s - x R ), with sample arm length x s and reference arm length x R : 

i D (k) oc I s (k) + I R (k) + ^/l s (k).I R (k) • 2cos(k • AOPL), (1) 

where k is the wavenumber (k = 2iz/X) and I s and I R are the signal intensities from the sample 
and reference arm, respectively. 

Since we use sdLCS for doing spectroscopy, we are interested in the spectrum Is(k). The 
sample arm contribution I s (k) to i D (k) is related to the sample's backscattered power spectrum 
S, which can be obtained from i D (k) using the procedure described in Section 2.2. Note that all 
bold faced characters in this article represent wavelength (k) dependent parameters. By 
adjusting x R , we can control the geometrical round-trip path length I in the sample at which 
we obtain S and by step- wise alteration of i, we obtain a data set S(i) [1-3]. Subsequently, we 
use a single exponential decay model to obtain the attenuation coefficient spectrum jii t of the 
sample: 

S(^) = ae M , (2) 

assuming local sample homogeneity for all path lengths in the data set S(i). The model is 
fitted to S(i) with free running fit parameters (a and |u t ), similar to our analysis for tdLCS [1- 
3]. The contributions of the scattering coefficient \i s and the absorption coefficient jn a can be 
obtained by fitting their individual contributions to the measured |u t , as described in Ref. [3]. 

The axial measurement range of the spectrograph is limited by its maximal measurable 
path length difference and sensitivity roll-off in depth [4]. Therefore, the spectra S(i) are 
acquired within a path length window Ai, which determines the spatial resolution for the |Li t - 
determination (see Section 2.3). Since a small axial measurement range is related to low 
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spectral resolution [11], a tradeoff exists between spectral and spatial resolution in sdLCS. 
However, we will show in this article that by selecting the right spectrograph, both 
sufficiently high spectral (6 nm) and spatial resolution (axial: 22 um) can be achieved. 




0 time Nt -1/(2t) 0 1/(2t) 0 time Nx 




Fig. 2. Signal acquisition and processing in sdLCS. See Section 2.2 for details. 
2.2. Removal of the DC and modulation terms 

The sample's backscattered spectrum S(i) in Eq. (2) cannot be obtained directly from i D (k), 
because i D (k) contains unwanted DC (i.e. non-modulated components) and modulation 
components (Eq. (1). Since the spectral features in S — originating from the sample's 
characteristic absorption and backscattering — may exhibit oscillation frequencies similar to 
the modulation on i D (k), S cannot be directly obtained from i D (k). Similar to the method 
described for OCT in Ref. [9], we will isolate S from i D (k) by phase modulation and 
frequency domain filtering of the signal. 

We start by modulating i D (k) using an oscillating mirror in the reference arm of our LCS 
system. This introduces a Doppler shift f D for every wavelength: f D = 2v R A, = 2AR-f R -k/7i, with 
v R the velocity, f R the oscillation frequency and AR the scanning range in one direction of the 
reference mirror. Acquisition of N spectra i D (k) over a time interval N(x + x d ) with integration 
time x and dead time x d results in a data set io(k,t), as illustrated in Fig. 2(a). Since x d (12 us) 
« x (6 ms) for the measurements described in this manuscript, we apply N(x + x d ) »Nx in the 
remaining part of our analysis. The frequency of the modulation on i D (k,t) changes as a 
function of time, because the movement of the reference mirror induces a change in AOPL 
(vertical direction in Fig. 2(a)). Figure 2(b) illustrates the time-modulation on i D (k,t) by f D for 
one wavelength (horizontal direction in Fig. 2(a)). Since the sensitivity of the measurement is 
largest around AOPL = 0 or 'zero delay' (see Eq. (3), the range of AOPL is chosen such that it 
crosses zero delay at /4Nx (Fig. 2(c)). Adequate sampling of the time-modulation on i D (k,t) is 
determined by the Nyquist criterion, i.e. the sampling frequency f s = 1/x needs to be larger 
than 2f Dmax (the largest value of f D within the investigated spectral range). For the integration 
time, this results in x < 7i/(2v R k max ) = A, min /(4v R ), with k max the largest wave number and X m \ n 
the shortest wavelength within the investigated spectral range. 

Fourier transformation on i D (k,t) with respect to time provides the frequency content of the 
signal i D (k,f), which contains the modulation frequency f D , its mirror image and a DC 
component (Fig. 2(d)). By filtering out the part of iD(k,f) containing only the positive 
frequencies + f D , the unwanted DC component and the negative frequencies -f D are lost. 
Inverse Fourier transformation on the filtered io(k,f) with respect to f provides a smooth, non- 
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modulated i D (k,t) that is confined in both amplitude and phase (Fig. 2(e)). Averaging of i D (k,t) 
over time results in the final spectrum S that we can use for the determination of \i t (Fig. 2(f)). 
A computational validation of this method will be shown by means of a simulation in Section 
3.3. 



2.3. The path length window Al in sdLCS 

In our approach for sdLCS, the path length window Ai over which S(X) is acquired results 
from a convolution of the reference mirror scanning window Ai R in the sample (Ai R = 2AR/n, 
with n the refractive index of the sample), with the path length window that is probed by the 
spectrograph AZ S . In practice, this path length window can be approached with Al = AZ R + 
Ai s (inset of Fig. 3). The spectrograph probing window Ai s depends on the maximal 
measurable path length difference by the spectrograph AOPL max (or imaging depth d max = 
V2AOPL max ) and the sensitivity roll-off with AOPL of the spectrograph. Due to the Nyquist 
criterion, AOPL max = 2n/8k with 'pixel resolution' 8k = (km^-kminVNp and N p the number of 
pixels of the spectrograph. The sensitivity roll-off of the spectrograph can be calculated using 
[11] 



I(AOPL) ■■ 



sin 2 (;rAOPL/[2AOPL max ]) 
(;rAOPL/[2AOPL m j) 2 



exp 



7i 2 (Ak/5k) 



81n2 



AOPL 
v AOPL , 

V max J 



(3) 



in which Ak is the spectral resolution of the spectrograph, assuming a Gaussian beam profile. 

For our sdLCS experiments, we use a spectrograph with a sensitivity roll-off that is two 
orders of magnitude smaller than AOPL max (Section 3.1). Therefore, we can approach Al s 
with the full width half max (FWHM) of the roll-off function: AZ S = FWHM roll . off /n. Hence, 
the path length window in sdLCS that determines the spatial resolution of the measurements is 
determined by: Al = Ai R + Al s = (2AR + FWHM roll . off )/n. 
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Fig. 3. Theoretical and measured sensitivity roll-off of the spectrograph. Inset: illustration of 
the path length window Ai in sdLCS (Section 2.3). 



2.4. Sensitivity/speed analysis 

The SNR of any shot-noise limited LCS system is given by 
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Here, £ is the detector quantum efficiency; hv is the photon energy and f B w is the 
electronic detection bandwidth of the system [5]. For a spectral domain (sdLCS) system, f BW 
= 111 [4]. For a time domain (tdLCS) system, f BW is determined by the imposed Doppler 
frequency f = 2v R A, = 2L scan>TO -f scan>TD IX of the scanning reference mirror [1], with v R its 
velocity, L scanTD = 2AR its scanning range in air and f sca n,TD = fR its scanning frequency. A 
commonly used comparison results in an improved SNR with a factor N p (number of 
spectrograph pixels) for sdLCS compared to tdLCS [4,5,12], provided that £, L scan and the f scan 
(i.e. acquisition speed) are equal for both methods. 

For a more generalized comparison between sdLCS and tdLCS, we take into account the 
differences in quantum efficiency, scan speed and scan range of the two methods (X 0 and X B w 
are equal since we employ the same light source): 



SNR SD = N p 



SD 



(2f R x) 



2AR 

AL maxy 



SNR TD , (5) 



with AL max the maximal measurable path length difference by the spectrograph after removal 
of the DC and modulation terms. 

3. Materials and methods 

In order to validate our data acquisition and analysis approach for sdLCS, we simulated an 
absorption measurement with the data acquisition parameters of our experimental setup. In 
addition, we measured attenuation spectra u t from 480 to 700 nm on a phantom with known 
absorption u a and scattering u s coefficient spectra and we separated the individual 
contributions of u a and u s from the measured u t . To compare sdLCS to tdLCS, we performed 
a signal to noise analysis for both systems and repeated the phantom measurement with the 
tdLCS system described in detail in Ref. [1]. 

3.1. System and acquisition 

The sdLCS system described in this article is identical to the tdLCS system with a low-pass 
filtered supercontinuum light source (SC430-4, Fianium Ltd., UK) described in Ref. [1-3], 
except for the detection end of the system, which consists of a spectrograph rather than a 
photodiode/lock-in amplifier combination. Hence, the multimode detection fiber is connected 
to the detecting spectrograph (USB4000, Ocean Optics, USA), with X^ n = 345 nm, A, max — 
1042 nm and N p = 3648. 

For both the sdLCS and tdLCS measurements, we controlled 1 (1 = 0- 2000 urn, with i = 
0 the sample surface) by translating the reference mirror in steps of 27 |im. By translating the 
sample in the axial direction, focus tracking of the spot size (r = 4.5 um) in the sample was 
achieved. At every £, back scattered power spectra S(i) were obtained over a path length 
window in the sample of Ai = 43 um for sdLCS (see Section 3.2) and Ai = 44 um for tdLCS 
[1]. For both sdLCS and tdLCS, S(i) was temporally averaged (N = 250 spectra per i) and 
corrected for the background (i.e. the LCS spectrum obtained from the non-scattering glass of 
the inner cuvette wall of the sample). The integration time was set at the minimal stable 
integration time of this spectrograph, x = 6 ms. To meet the sampling frequency requirement 
(Section 2.2), the velocity of the reference mirror was set at v R = 12 |im/s (frequency f R = 0.3 
Hz, scanning range AR = 20 urn). The acquisition settings for the tdLCS measurements were 
identical to the settings in Ref. [1-3]. For both sdLCS and tdLCS, the most important 
acquisition parameters have been summarized in Table 1 . Note that the acquisition speed for 
one temporally averaged spectrum S(i) is 3.6 times faster for sdLCS. 

As described in Section 2, fitting the single exponential decay model S(i) = a-Qxp(-\i t -i) 
(free running fit parameters a and \i t ) to the sdLCS and tdLCS acquired S(i) vs. i, results in a 
\i t spectrum for both detection methods [1-3]. Uncertainties in a and \i t are estimated by their 
95% confidence intervals (c.i.). Since the spectral resolution in sdLCS (AX = 6 nm) is higher 
than the pixel width of 0.2 nm, S(i) was binned into wavelength regions of 6 nm for the 
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sdLCS measurements. The values of jn a and |Li s were obtained by fitting their individual 
contributions to the measured jii t , as described in Ref. [3]. 



Table 1. Acquisition settings for sdLCS and tdLCS 



Acquisition parameter 


sdLCS 


tdLCS 


Ai 


43 urn 


44 um 


AX 


6 nm (all X) 


4 nm @ X = 480 nm 
9 nm @ X = 700 nm 


Vr 


0.006 mm/s 


1.84 mm/s 


f R 


0.3 Hz 


23 Hz 


AR 


20 um 


40 um 


N = # averages per i 


250 


250 


Acquisition time of S(Z) 


Nx=1.5s 


N/(2f R ) = 5.4 s 



3.2. System characterization 

Using the specifications of the spectrograph and the theory described in Section 2.3, we 
calculated that AOPL max = 1880 um. After measuring the spectral resolution of AX = 6 nm 
with a 543 nm and 633 nm HeNe laser line, we used Eq. (3) to calculate that the FWHM roll . off 
is 18 jim ('theoretical' roll-off in Fig. 3). This FWHM roll . off corresponds well with the 
experimentally measured roll-off function, measured from the reflection of a glass slide in the 
sample arm as a function of AOPL between the sample arm and the reference arm (Fig. 3). 
Since the FWHM roll . off of the spectrograph is two orders of magnitude smaller than AOPL max , 
the path length window that is probed by the spectrograph can be approximated with Ai s = 
FWHM roll . off /n =13 urn. Consequently, the total path length window that is probed in our 
sdLCS measurements is Ai = Ai R + Ai s = (2AR + FWHM roll . off )/n = 43 um, with a fixed AR 
of20um. 

Based on Eq. (5), we can calculate the expected difference in SNR between both systems. 
The sdLCS spectrum covers approximately N p = 1150 pixels on the USB4000. Based on a 
maximum number of 60 photons/count and a pixel well depth of 1.1 0 5 electrons given by the 
manufacturer, the quantum efficiency (QE = #electrons/photon) of the 16-bits USB4000 is 
2.5% at 600 nm. Using the photodiode response of 0.43 AAV (2001, New Focus, USA), we 
can calculate that our tdLCS detector has a QE of 89% at 600 nm. With x = 6 ms for our 
sdLCS system, the sensitivity advantage of the sdLCS measurements with the current 
equipment reduces to 0.38. Hence, when using the USB4000 as a detector for sdLCS with the 
settings in Table 1, we can expect a decrease in SNR of 4.2 dB at 600 nm with respect to 
tdLCS. 

In order to compare the sensitivity of the current sdLCS system to the tdLCS system, we 
measured the SNR of the two detection methods, using a mirror behind an attenuating (2.0 
optical density) filter in the sample arm. The optical power reference arm was optimized for 
(near) shot noise limited detection and the acquisition settings in Table 1 were applied. For 
both sdLCS and tdLCS, the spectrally resolved SNR was determined from the measured 
spectrum S at AOPL = 0 using [5]: SNR = 10- log(S/var [noise]). Here, var[noise] denotes the 
variance of the noise spectrum, measured with a AOPL large enough to observe no 
interference signal within the path length window Ai. For both detection methods, the 
measured SNR is around 75 dB per wavelength bin between 480 and 550 nm. After 550 nm, 
the SNR of tdLCS increases to a value of 85 dB at 700 nm, which can be explained by the 
larger quantum efficiency of the photodiode in tdLCS for longer wavelengths. The SNR of 
sdLCS remains spectrally flat over the full wavelength region of 480 to 700 nm, resulting in 
an SNR of -3 dB lower compared to tdLCS at 600 nm, which is close to the predicted value 
of4.2dB. 

3.3. Simulation 

To validate the analysis algorithm described in Section 2.2, we simulated an absorption 
measurement on a sample containing oxyhemoglobin [13]. With this simulation, we tried to 
mimic the reality of the acquisition as closely as possible; therefore we incorporated 
summation over the integration time for every spectrum while simulating reference mirror 
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movement, a discrete wavelength axis and a limited spectral resolution. As a consequence, 
dynamic effects on such as fringe washout are taken into account. The contributions of 
scattering and noise were neglected in the simulation. 

The detecting spectrograph was simulated with the properties of the spectrograph used for 
our measurements, as well as the values for the integration time and reference mirror velocity. 
The source spectrum S 0 was simulated as a Gaussian with center wavelength Xq = 550 nm and 
bandwidth A, FWH m =140 nm. To account for the spectral resolution (AX) of the spectrograph, 
we convolved the spectrum with a Gaussian with a full- width at half maximum of 6 nm. The 
simulated sample arm power fraction was r| s = 0.1, resulting in a reference arm power fraction 
of r| R = 0.9. Within each integration time, 16 time-separated spectra at different reference 
mirror positions were added to simulate an acquired spectrum. A total number of N = 125 
acquired spectra was simulated over a time interval of Nx = 0.75 s. The initial AOPL was set 
at -4.5 jum to achieve zero delay at /4Nx. 

Since scattering was neglected in the simulation, the attenuation coefficient \i t of the 
sample can be simplified to the absorption coefficient |u a . The spectrum S(l) was calculated at 
i x = 2 mm and i 2 = 4 mm inside the sample, corresponding to depths di = 1 mm and at d 2 = 2 
mm, respectively. We obtained the absorption coefficient of the sample using Beer's law (Eq. 
(2): Ma = -[V 2 (i 2 - i^T'MSdi)/ S(i 2 )l 

3.4. Phantom 

We prepared a phantom, consisting of 0.096 vol% NIST-certified polystyrene spheres (0 602 
+ 6 nm, Thermo Scientific, USA) and 10% magenta dye (Ecoline #337, Royal Talens, The 
Netherlands). The |n s of the polystyrene spheres was calculated using Mie theory and 
integrated over the size distribution of the spheres (2*SD), as described in Ref. [2]. and the jn a 
of the dye was determined in a separate transmission measurement on the dye only, as 
described in Ref. [1]. The measured tdLCS spectra are broadened by Doppler frequency shifts 
induced by the Brownian motion of the polystyrene spheres. Therefore, both reference spectra 
for jli s and jn a were convolved with a Lorentzian with a line width of 9 nm for adequate 
comparison to the LCS spectra [1-3]. For the sdLCS measurements, Brownian motion does 
not influence the spectrum S, because the Doppler shifted frequencies remain within the 
bounds of the frequency domain filter (Fig. 2(d)). Since the reference spectrum of |Li a was 
obtained with the same spectral resolution as the LCS spectra, we did not broaden this 
spectrum. The reference spectrum of jli s was convolved with a Lorentzian with a line width of 
6 nm (AX of the spectrograph). 

4. Results 

4.1. Simulation 

Figure 4 shows the result of the simulation described in Section 3.3. The i D (X) at i\ and i 2 
show little influence of the sample's absorption (Fig. 4(a)). However, after frequency domain 
filtering of in(X), the resulting spectra S at i x and i 2 clearly show the presence of the 
oxyhemoglobin absorption peaks (Fig. 4(b)). From Fig. 4(c), we can conclude that our method 
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Fig. 4. Simulation of a |Li a measurement in sdLCS. a.) input spectra i D and b.) filtered spectra (S) 
at path lengths ii and i 2 inside the sample, c.) input |Li a and recovered |n a . 
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of data acquisition and analysis for sdLCS fully recovers the input \i a of this simulation. Even 
though the spectral features of the oxyhemoglobin absorption exhibit oscillations with a 
frequency comparable to that of the modulation on i D (A), |u a is recovered without any artifacts. 

4.2. Phantom 

Figure 5 shows the measured |u t spectra by sdLCS and tdLCS on the polystyrene-dye 
phantom. The \i t spectra from both detection methods agree within the estimated uncertainties 
(error bars) for nearly all wavelengths. Also the fits on \x t and the |Li s -contributions to the fits 
agree well, and the |n s -contributions are in good agreement with the Mie calculated \i s . The 
fitted dye concentrations were 10.1 + 0.4% for the sdLCS measurement and 9.5 ± 0.4% for the 
tdLCS measurement, which are both very close to the expected dye concentration of 10%. 
Note that part of the difference in spectral shape is caused by the presence of Brownian 
motion influences for tdLCS, versus the absence of it for sdLCS (Section 3.4). 

15 1 sdLCS 
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Fig. 5. Comparison of an sdLCS measurement to a tdLCS measurement of jn t , Ha and jli s on a 
polystyrene-dye phantom. 

5. Discussion and conclusion 

In this article, we demonstrated and validated a new detection method and analysis algorithm 
for spectral domain detection in LCS. We showed that this method can be used for measuring 
local u t , u a and u s spectra in turbid media. The results of the sdLCS phantom measurement 
agreed well with the expected optical properties of the phantom and were comparable to the 
measured optical properties by tdLCS (Fig. 5). Whereas we did not explicitly show the 
localization aspect of our method by measuring on homogeneous media, previous studies with 
tdLCS proved that localized measurements of optical property spectra are possible within 
distinct tissue volumes [3]. 

The reason to validate this new detection method, is its theoretical sensitivity advantage 
compared to time domain detection (Section 2.4). As a consequence, the acquisition speed of 
the measurement can be reduced with respect to time domain detection, without reducing the 
signal to noise ratio. The sdLCS system that was validated in this article was more than three 
times faster than the tdLCS system, but not more sensitive (Section 3.2). A sensitivity 
advantage for sdLCS could not be demonstrated, because the quantum efficiency of the 
spectrograph is not optimized for these experiments. However, if the quantum efficiency 
would be similar to that of our tdLCS detector, the sensitivity advantage would be 1 1 dB at 
equal acquisition times for tdLCS and sdLCS (Eq. (5). 

An additional shortcoming of the used spectrograph is its relatively slow acquisition speed 
(167 spectra/sec). Slow acquisition speeds result in a decrease of the modulation depth in 
i D (A) ('fringe wash-out'), due to the reference mirror movement induced change of AOPL 
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within the integration time. As a consequence, an additional decrease in sensitivity can occur, 
accounting for approximately 0.5 dB in our sdLCS system. Hence, not only higher quantum 
efficiency, but also faster acquisition speed is required to further improve the performance of 
our sdLCS system. Since high quantum efficiency line scan cameras are available with line 
rates up to 140 kHz, enhanced sensitivity and/or acquisition speed are realizable for sdLCS. 

As we showed before for tdLCS, it may be advantageous to support the selection of the 
volume for optical property determination with an image of the investigated tissue volume [3]. 
Similar to tdLCS, this image can be reconstructed from the backscattered intensity as a 
function of depth. Our method for removal of the DC and modulation components in sdLCS 
also removes the 'mirror image' (complex ambiguity) when AOPL = 0 is situated inside the 
sample. As a consequence, this method for sdLCS can be used to reconstruct an image 
without cross-talk between i D (AOPL) and i D (-AOPL) within the investigated path length 
window inside the sample. 

The only technique that is reported to have comparable performance to our sdLCS system 
in terms of localized measurements of optical properties, is dual window sOCT that was 
developed by Robles et al. [8,14]. Whereas Robles et al. report that their method achieves 
higher resolution in both the spectral and the spatial domain by avoiding the inherent tradeoff 
between the two, their system may suffer from the effects of the sensitivity roll-off of the 
detecting spectrograph and the absence of focus tracking, which influence measurement 
depth. An advantage of our sdLCS system is that it does not rely on refractive index- 
dependent corrections for the unwanted signal attenuation due to these effects. This facilitates 
the exact determination of jli s and |u a contributions to the measured \i t . Since the current spatial 
and spectral resolution of our sdLCS system are comparable to the tdLCS system (Table 1), 
the resolution of our system is sufficient to obtain relevant determinations of local tissue 
chromophore concentrations, e.g. within the dermis, or epidermis only as shown in Ref. [3]. 

In conclusion, we have demonstrated and validated a new approach for sdLCS, which has 
high potential to improve the accuracy and speed of localized optical property measurements 
by LCS. Undoubtedly, this will lead to improved clinical utility of the technique, e.g. for the 
non-invasive determination of blood composition (hemoglobin/bilirubin concentration and 
oxygen saturation). 
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